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Abstract 

Level curvature is a measure of sensitivity of energy levels of a 
disordered/chaotic system to perturbations. In the bulk of the spec- 
trum Random Matrix Theory predicts the probability distributions 
of level curvatures to be given by Zakrzewski-Delande expressions [F. 
von Oppen Phys. Rev. Lett. 73 798 (1994) & Phys.Rev. E 51 2647 
(1995); Y.V. Fyodorov and H.-J. Sommers Z.Phys.B 99 123 (1995)]. 
Motivated by growing interest in statistics of extreme (maximal or 
^5 ' minimal) eigenvalues of disordered systems of various nature, it is 

natural to ask about the associated level curvatures. I show how cal- 
culating the distribution for the curvatures of extreme eigenvalues in 
GUE ensemble can be reduced to studying asymptotic behaviour of 
orthogonal polynomials appearing in the recent work C. Nadal and S. 
N. Majumdar J. Stat. Mech. 2011 P04001 (2011). The correspond- 
ing asymptotic analysis being yet outstanding, I instead will discuss 
solution of a related, but somewhat simpler problem of calculating the 
level curvature distribution averaged over all the levels in a spectral 
window close to the edge of the semicircle. The method is based on 
asymptotic analysis of kernels associated with Hermite polynomials 
and their Cauchy transforms, and is straightforwardly extendable to 
any rotationally-invariant ensemble of random matrices. 
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*The text is based on the presentation at the 5th Workshop on Quantum Chaos and 
Localization Phenomena, May 20-22, 2011, Warsaw, Poland 



1 



1 Introduction 



Let if at stands for N x N random Hermitian matrix belonging to a certain 
invariant ensemble which is characterized by the joint probability density of 
N real eigenvalues Aj, , i — 1, . . . , N of the form 



i ^ 
7MA 1; A 2 , . . . , X N ) = i- e-f^ TT (A n - A,) 2 . (1) 

in terms of the potential ^(A), with being the appropriate normaliza- 
tion constant. In particular, the simplest choice V(\) = A 2 corresponds 
to the so-called Gaussian Unitary Ensemble (GUE) whose mean eigenvalue 
density is given in the limit N — > oo by the Wigner semicircle law p(/i) = 

jf (zuLi - K)) = ^:a/4 - // 2 , |A| < 2. This law shows, in par- 
ticular, that typically the minimal X m in and the maximal \ ma x eigenvalues 
approach ±2, respectively. For large but finite N one of the most important 
chracteristics of GUE spectrum appears to be the Tracy- Widom (TW) law 
for the distribution of this extreme eigenvalues given by pQ 



V(\ max ) = F 2 ( A "^ 2/3 2 ) , F 2 (x) = ^ exp 
where q(z) satisfies the Painleve II equation: 



(z — x)q 2 (z) dz 



(2) 



q"(z)=2 q Z(z) + zq(z), q(z ^ oo) ~ Ai{z) ~ — =r TJl e"^ (3) 

Similar distributions are also known for other symmetry classes of random 
matrices. TW distributions are highly universal, retain their validity not 
only for invariant ensembles with a quite general potential V(X), but also for 
a very broad class of Hermitian random matrices with independent entries 
[2], and emerge in several disordered/chaotic physical systems. E.g they are 
relevant for describing directed polymers in disordered media [3], spectral gap 
fluctuations in disordered metal grains in contact to a bulk superconductor 
jl], fluctuations of output power in coupled fiber lasers [5], etc. 

Consider now a general perturbation + jW, where 7 is the control 
parameter, and W is a fixed given matrix. One can pose a natural ques- 
tion of characterizing sensitivity of the minimal/maximal eigenvalue A m j n (7) 
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to a generic perturbation by considering the standard perturbation theory: 
(7) = ^min + 7 V + 7 2 C + . . ., where 



Here |n), A n for n = 1, . . . , iV is the set of of eigenvectors/eigenvalues of H^, 
that is H\n) = X n \n), and |m) stands for the eigenvector corresponding to 
the unperturbed minimal eigenvalue A min . The coefficients V and C in (j4j) 
are frequently called in the physical literature the "level velocity" and the 
"level curvature", respectively. This terminology is inherited from the use 
of eigenvalues of random matrices as a model of highly excited energy lev- 
els of disordered/chaotic quantum systems, see e.g. [6\. As is well-known, 
the components of eigenvectors of invariant random matrices are statisti- 
cally independent from the eigenvalues, and also in the large- N limit behave 
essentially as independent, identically distributed Gaussian variables with 
variance 1/N. This makes calculating the distribution of V within the Ran- 
dom Matrix Theory (RMT) context a straightforward task. At the same time 
finding statistics of the level curvature C is a much less trivial problem. Var- 
ious aspects of the level curvatures for eigenvalues in the bulk of the spectrum 
of random matrix ensembles , as well as for disordered and chaotic systems 
attracted quite a considerable interest in mid- '90s. The corresponding cur- 
vature distributions were conjectured by Zakrzewski-Delandep] on the basis 
of numerical simulations, and derived in the limit N ^> 1 for Gaussian en- 
sembles in [8] and independently by a different method in [9] (see yet another 
technique in [10]). One expects the results to be universal, that is to hold 
for a broad class of random matrices sharing the same global symmetries. 

To this end it seems natural to pose questions related to sensitivity of 
extreme eigenvalues to perturbations. Apart from a generic interest, from a 
somewhat different angle the above expressions characterize sensitivity of the 
ground state of the so-called spherical spin-glass model[TT] to perturbations 
in random interactions. In addition, one can show also that in the framework 
of the same spin-glass-type spherical model the so-called nonlinear suscep- 
tibility of the ground state with respect to external magnetic field can be 
reduced to a similar, but a more complicated expression. Those observations 
provide an additional motivation to try to develop regular tools for statisti- 
cal characterization of level curvatures for extreme eigenvalues. Here I report 
some preliminary steps in this direction. 
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The presentation below will have the following structure. I will start with 
recapitulating a simpler problem of the level curvature distribution, with the 
curvatures being sampled over all levels in a small spectral window around 
some point in the spectrum. Departing from the specific methods used in 
earlier papers j8], j9] for Gaussian RMT ensembles, I show that for any rota- 
tionally invariant ensemble of Hermitian random matrices the problem can 
be reduced to the asymptotic analysis of the orthogonal polynomials, and 
their so-called Cauchy transforms. After that I will go back to the problem 
of sensitivity of extreme eigenvalues and show how the corresponding curva- 
ture distribution can be formulated in terms of a special class of orthogonal 
polynomials studied recently in [T2]. The full asymptotic analysis of the re- 
sulting kernel is still outstanding, I will briefly show instead how the new 
formulation allows to reproduce known GUE results in the bulk, and then 
concentrate on the GUE "soft-edge" scaling regime which was not yet studied 
before. Appendices A,B contain some technical details related to asymptotic 
analysis of Hermite polynomials. 

2 Level curvature distribution for invariant 
ensembles: general consideration 

2.1 Spectrally-averaged curvature distribution in terms 
of orthogonal polynomials 

The probability density of level curvatures averaged over all the eigenvalues 
around the point \x in the spectrum is defined as 



where (...)# stands for the ensemble average, p(fi) = (J2 n =i ~ ) 

\ / H N 

stands for the mean eigenvalue density around the point //, the curvature 
related to the m-th eigenvalue being defined as 
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and C typ stands for the typical curvature scale identified via considering the 
typical contribution of two neighbouring eigenvalues 



C ty P = ^(m|W|7i)(7i|W|m) = irp(p) y typ , y typ = ^TrW 2 (7) 

with the parameter A = l/(Nnp(p)) defining the mean spacing between the 
neighbouring eigenvalues. The bar stands for the average over eigenvectors, 
and we assumed above that both \n) and \m) are independent iV— component 
gaussian complex- valued vectors with mean unit length (see below). In what 
follows we also assume that a generic perturbations is such that the vari- 
able ytyp = 0(1) for N — > oo. More precisely, we consider only full-rank 
perturbation W such that all its N eigenvalues are of the order of unity. 

Let us note for the future use that from the definitions (I5ll6l) the mean 
curvature for levels around a point /i in the spectrum is given by 

i r°° d\ 
(c m ) H = j^^—x (8) 

where we introduced the standard eigenvalue two-point cluster function 

/>oo 

n 2 (ji,\) = N(N-l) / V N (X 1 = fi,X 2 = X,X 3 ,...,X N ) d\ 3 ...X N (9) 



Note that the eigenvalues are strongly correlated only over the distance com- 
parable with typical level spacing A which is negligible in comparison with 
the total length of the spectral support (assumed to be for simplicity a sin- 
gle interval). It is then easy to see that the leading order result is obtained 
by neglecting the correlations and using 7£ 2 (/i, A) ~ N 2 p(p)p(X) in (jSj). In 
particular, for GUE the mean curvature in the limit N — > oo is simply given 
by 

f 1, M<2 

(Cm) H = Vtyp < (10) 
I 2 , |//| > I 



where we have assumed that /i is fixed when iV — > oo, and used the known 
formula for the averaged resolvent of the GUE matrix. We also included for 
further reference second line corresponding formally to the situation when 
the observation window is chosen to be outside the support. 

After this digression let us proceed to calculating the Fourier transform 
K,{uS) = J V(c, fi)e~ luJC dc of the probability density for normalized curvatures 
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c = C '/C t y P . Introducing for notational shorthand the vector \W m ) = W\m), 
we have after straightforward manipulations: 



u (n\W m )(W m \n) 



m=l n(^m) ^ * OT , ,, 

(11) 

We will perform the ensemble average in steps, and start with averaging 
over the eigenvectors \n) with n — 1, 2, . . . , m — 1, m + 1, . . . , N. In doing 
this we assume that in the large— N limit different eigenvectors are effec- 
tively statistically independent for different n, and moreover the components 
Hi = (i\n) of a given eigenvector in any basis \i) can be treated in the same 
limit simply as i.i.d. complex Gaussian-distributed numbers with variance 
1/N. Actually, one can relatively easily verify that if one takes any finite 
subset of I eigenvectors such that I is fixed when iV — > oo then the required 
properties follow. The latter statement is rigorously proved in the mathemat- 
ical literature, see e.g. [13]. Therefore, our method implicitly assumes that the 
number I of terms which effectively contribute to the curvature defined in 
is much smaller than N. This is very plausible in view of the denominators 
growing roughly linearly with n, but strictly speaking remains a conjecture. 
Using such an assumption we can easily perform the eigenvector average for 
a given \n), which simply amounts to using the Gaussian identity: 

— I e - N ( n \ n )- ia "( n \ w ™)( w ™\ n ) d 2 \ n } ' 



M J 1 ' det(l + if \W m )(W m \) 

1 « n = l — (12) 



l + i^(W m \W m y Cty P fi-X n 
This immediately allows us to write 



N N 



An 

where the brackets now stand for the averaging over the joint probability 
density (CQ), and we have introduced (random) variable y m = ^ Wr ^J^ = 
(m\w |m) ^ anc j exploited the Hermiticity W = W^), with the bar standing for 
the remaining average over that variable. 
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Let us briefly discuss how one performs the latter average by calculating 
the distribution of the variable y m . Exploiting the mentioned Gaussianity of 
the individual eigenvector we immediately find for the corresponding char- 
acteristic function: 



oil iw^— Trw 2 9_Tny 



exp { ^ W r W f = det(i-^ w ») ^ e " ' 



(14) 

As for a "generic" perturbation we must have TiW 2p = O(N) for all inte- 
ger p > 1 one can for large N effectively retain only the first term in the 
exponential in ffTTj) . We conclude that the distribution of y m is effectively 
S— functional in the limit iV ^> 1, so one can replace y m with its average 
value: y m ->y^ = TrW 2 = ^ = AN. This gives 



JV-l 



]C(u) = -^{5(jjl- X N ) ff ~, Xn a > (15) 

where we have used that by permutation symmetry of the joint probability 
density (JTJ all the iV terms in ( 1T31) produce identical contribution upon av- 
eraging, so it is sufficient to consider averaging of a single term with m = N. 
To perform the remaining averaging over eigenvalues it is convenient to in- 
troduce the ensemble of (N — 1) x (N — 1) Hermitian matrices H^-i whose 
joint probability density of eigenvalues is given by (cf.(JTJ ) 



1 N ~ l 
V N -i(Xi, A 2 , . . . , A^_0 = ^— e-Wi' v ^ IT ( A « - A ') 2 • ( 16 ) 



n<l 

This allows us to rewrite (JTJ as 

N-l 



V N (Xi, A 2 , . . . , Ajv) = -^ e - NV ^ TT (Ajv - A n ) 2 ^-i(Ai, A 2 , . . . , A^-i) 

ZjV n=l 

(17) 

Introducing the characteristic polynomial det (/i — i/jv-i) one can notice that 
the integration over Ai, . . . , Aat_i in (TT5|) amounts to evaluating the following 
object 

'[det^-iWV _ 

|('V-ll" 

H N - 



det (e - Fjv-i) / [c^v-i] 2 ' 
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where we have denoted e = \i + iu> A, A = wp ^ N and the constant co- 
efficients [cat_i] 2 will be defined below. In terms of the above function the 
Fourier-transform of the curvature distribution is simply given by 

J V(c, f x)e~^dc=^^e-^ 2 j^^^ + l u J A) (19) 

The correlation functions of product/ratios of characteristic polynomials of 
random matrices, like that featuring in ffl8l) . were considered in full gener- 
ality in [HI [15]. There it was found how to express these objects in terms 
of orthogonal polynomials generated by the potential V(x). Namely, intro- 
ducing a measure on the real line as djj,{x) = e~ NV<yX ^ dx, one can define the 
unique set ^(x) of associated monic orthogonal polynomials satisfying 

n k (x)nj(x)dfi(x) = [c k ] 2 S jk (20) 

As is well-known (see e.g. [To 7 } [T7] and references therein) the product of 
the coefficients [c^] 2 gives the normalization constant in ([1]): Zjv_i = (N — 
1)! rifc^o 2 l c kf- Further, define the so called Cauchy transforms 

hk{e) = ^- ! ^ dfi(x), Ime + (21) 
2m J x — e 



In particular, in [15] an expression was derived relating general averages like 
those featuring in (IT81) to the determinants of 2 x 2 matrices whose entries are 
certain bi-linear combinations of the polynomials and their derivatives (the 
so-called "kernels"). After specifying the general formulae for our particular 
case the correlation function defined in (fl9|) takes the following explicit form 



W!{fi,fi) - -(e - fJL)W 2 (n, fi) 



+ (e-n)F 2 (n,e)W 1 {n,n), 

(22) 

in terms of the following kernels: 

Fi(/i, e) = h N (e)'K N ^ l (p) - h N -i{e)n N (n) (23) 
F 2 (fi, e) = hNieWN-iit*) - h N ^i(e)n' N (fi) (24) 

and 

Wti/i,/!) = ii' N (fi)n N ^i(fi) - n N (fi)n' N ^ 1 (fi) 
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So formally the problem amounts to finding asymptotic approximations 
for the orthogonal polynomials and Cauchy transforms for a given potential in 
the specified spectral regime. Various techniques are available for performing 
such an analysis, the Riemann-Hilbert approach (see e.g. [16] ) being the most 
powerful, especially for proving the universality of the required asymptotics 
for a broad class of rotationally-invariant ensembles. We will not pursue 
this line here but rather show later on how the above formulae reproduce the 
known Zakrzewski-Delande expressions in the bulk of the spectrum < 2 of 
the Gaussian Unitary ensemble in the large- N limit, and then study the soft- 
edge case. Before doing that we however come back to addressing the original 
problem of the curvature distribution for extreme eigenvalues, and formulate 
it in terms of asymptotics of a special class of orthogonal polynomials. 

2.2 Curvature distribution for extreme eigenvalues: or- 
thogonal polynomial formulation 

Consider again the perturbation of the extreme eigenvalue {\ m in for definite- 
ness) with the curvature given by (jlj) and normalized to the typical curvature 
at the soft edge C typ = N~ l ^y typ (see Sec. 3.2 below): 

/ ' ' (m|jy \n) (n\W\m) 



\m=l \ L HP n^m 




N( S[c 



1 ^ (l W n (n W l ) \" \ 
7T- \ \ 1 1 X(a„>ao ) (27) 

C typ n =2 Al " A ™ J n=2 / Hn 

where averaging goes over the joint probability density ([T]), and we have in- 
troduced the indicator function: X(A) = 1 if A is true and zero otherwise, and 
exploited the permutation symmetry of ([1]). Introducing the corresponding 
Fourier-transform (known as the characteristic function) K, m {oj) = (e tu)C ) and 
averaging it over the Gaussian eigenvectors |n), with n = 2,3, ...,N, and 
then over the remaining eigenvector |1) corresponding to \ m in yields, in full 



analogy to ([T5 



N \ - \ N \ 

^)=(n Al _r + u n^>M . ^ 
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At the next step we introduce an ensemble of (N — 1) x (N — 1) random 
Hermitian matrix Mjv-i with the eigenvalues X 2 , . . . ,Xn and the (normalized) 
measure 

N N-l 

Px min (M N -i) = „ e-T^^Hx^) H(K-\ P ) 2 (29) 

v ' n=2 n<p 

where Zjv-i(Ai) is the appropriate normalization constant. This distribution 
(and its normalisation) depends on Ai(= \ m in) as & n external parameter. In 
terms of such an ensemble we easily see (cf. ffTSl) ) 

M«) = ^ r d\ min e-^ min ) ZN _ i{Xm . n) I TT !^^L\ 

«/— oo \n=l min ~ n / M 

(30)" 1 

where the averaging (...)_j if goes over the probability density (1291 . 

We conclude that the analysis of the above expression amounts to study- 
ing the orthogonal polynomials generated by the measure P\ min {M^^i). In- 
deed, we can introduce the measure dn\ min (x) = e~~ v ^X(x>\ mm ) dx, with 
^k(x] A m j n ) standing for monic orthogonal polynomials with respect to this 
measure satisfying 

/nhW^U*uH = MW** (3D 

The above polynomials depend on A m j n as external parameter, and in this 
way the analysis of the level curvature distribution for the minimal eigenvalue 
amounts to extracting the large-N asymptotic behaviour of such polynomials 
and their Cauchy transforms for N — > oo. Several important steps in such 
an analysis for the Gaussian case V(x) = x 2 were reported recently by Nadal 
and Majumdar in [12J, but further work is needed to include the Cauchy 
transforms into considerration to be able to extract the ensuing curvature 
distribution in the explicit form. The problem is non-trivial and is currently 
under investigation [T5]. Below we return to considering a somewhat simpler 
case of the spectral-averaged curvature distribution, both at the edge and in 
the bulk. 
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3 Spectral-averaged GUE curvature distribu- 
tion 



3.1 Bulk of the GUE spectrum 

Denote Tik{x) = Pk(x) = x h + . . . the monic orthogonal polynomials w.r.t. 



AT 2 



the standard Gaussian measure on the full line djjx min {%) = e 2 x dx, that 
is 

J p k (x) Pj (x)e-^ x ' 2 dx = c 2 k 5 jk , c 2 k = j^\J^ (32) 

Those are actually the classical Hermite polynomials. Renaming the spectral 
parameter fj, — > x, and y = so that A = our goal is to calculate (see 



/N ~ 1 u 

^'^'W^WW^f^ Hx ' x + 'n v) (33) 

Relevant bulk asymptotic expressions for Hermite polynomials, Cauchy trans- 
forms, and the kernels involved in the curvature distributions are well-known, 
but to make the present text self-contined I recover them in the Appendix A 
directly from the integral representations. Parametrizing a bulk point of the 
spectrum x = 2 cos</> G (—2, 2) for i], ( of the order of unity and iV > 1 we 
have 

27Ti TP ( , V , C 



c N-l 



Fl ( x + l jX + j-\^ e -f (t-v)+iMx)(t-v>c (34) 



where we denoted = sign/m(£). 



Taking into account — ^-4- — y ie N we conclude that the required large— N 



2ni 



asymptotics for any real u> and y > of the order of unity is given by 

—11} m i P -N p -iu%y-Kp{x)\w\y 



F 1 (x, x + 1— yj » i e- w e- w S v-'VWMtf (35) 

whereas in the view of the exact relation F 2 (x, x + ijjy) = N-^ ^=0-^1 (x + jf,x 
we have 

F 2 (x,x + i^y) » N |-is w vrp(x) Fx (x, x + i^-yj (36) 
We also have 

Wi(x,x) « 2A^e~ iV 7rp(x)e^ 2 , ^(1,^11^(1,1) (37) 
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Substituting all this to (122]) with p = x, e = x + z^y we find that 



-N 2 / . W 

e 2 y- ^r, x + 



1 . w 

1 1— y J\ x 

2 N y 



■ u 



x 



- 1 s^n 



p(x) 



= 2i Ne~ 2N irp(x) e~ iuj % v-M*)Mv(i + y\uj\np(x)) (38) 

Now we substitute all this into (133]) and use the Stirling approximation 
(jV - 1)! = ^^N N e~ N . Finally we arrive at 



(39) 



This is indeed exactly the expression leading after the Fourier-transform to 
the Zakrzewski-Delande curvature distribution for GUE ensemble, see [TJ El 



2 oc x 

p( c x) = -2 , k = 7rp(x)y = 1, c = - y = - — n 

71 [(c - c ) 2 + k \ 2 2vrp(x) 



(40) 



Note: The above calculation can be further shortened if we first notice 
two useful identities, the first one being 



N P (p) = ^e~^ 2 ([det (p - Mjr-x)] 2 )^ 



(41) 



and second one Wi(p,p) = ([det (/i — M/v-i)] 2 ) M ^ (see e.g. [H]). When 
combined together they produce the following relation 



-e 2^ Wi (//,//) = 1 



Np(p) Z 



N 



(42) 



which after being substited to (122"]) helps to rewrite the Fourier-transformed 
curvature distribution in the most concise form: 



V(c, p)e lu)C dc = 7 [p,p + i—y 



(43) 



12 



where we have defined 
J"(/i,e) = Fi(/i,e) 



2 [e ^VO^J 



+ (e-^)F 2 (^,e) 



(44) 



and Fi )2 (/i, e) 



2~; 



^i, 2 (/i,e). 



In particular, such a form turns out to be more convenient for extending 
the calculation to the spectral-averaged curvature distrubution at the "edge 
of spectrum" regime which is to be considered in the next section. 



3.2 Soft edge of the GUE spectrum 



Consider the "soft edge" regime p = x = 2 + where ( G (—00, 00). The 
mean density of eigenvalues p(x) in this regime is well-known and scales as 

1 



p{x) 



JVV3 



p((), with p(0 = ^'(C) 2 -^(C)^"(C) 



(45) 



where Ai(() stands for the Airy function. The corresponding mean level 



spacing is then A 



where y 



Relevant soft edge asymptotics 



of Hermite polynomials, Cauchy transforms, and the kernels involved were 
considered, for example, in [19] (see also [20], [21] and [H]). For completeness, 
we reproduce them in detail in the Appendix B, see in particular expressions 
(1T591 for and also e q uations (HZU), (HZOD for F 2 (x, x + j^) and 

Fx (x, x + 



2Wi(x,x) 

j^jw)- Remembering e—p — i 

c 



algebra to 

V(c,p)e- wc dc = T\ 2 



-^273 we reduce (I43II44P after a simple 
( + iwy\ ,_„_iivi/3 



AT2/3' 



2 + 



AT2/3 



oc — me 



*(C,w)-i*(C,w)^lnp(C) 



where 



(46) 
(47) 



$(C, w) = a((, u>) Ai'(C) - «'(C, w) Ai(C), *(C, w) 
= «(C, w) Ai (C) - 01 (C, w) Ai (C) with a; (C,a;) = — a(C,a;), etc. 

The exponential factor e ~ lNl/3u) y simply fixes the constant shift in the 
curvatures c = C/C typ to be cq = iV 1 / 3 ?/ = (remembering the correspon- 



dence y = iV 1 / 3 ?/ this value coincides with the mean bulk value Cq 



111 
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the edge limit x — > 2). To get rid of such a shift and of the related extra p(Q 
dependence in typical curvature we redefine "shifted and scaled" curvatures 
c sc = = 7rp(C)( c — c o) for eigenvalues around the point x = 2 + ^573. 
This will allow us to omit the factor e ~ lNl/iui y an d set y — 1 in the definition 
(TWDfor a(C,o;). 

Note that by setting in (jTUJ) /x = 2 + ^§73 yields for the averaged curvature 

(C) /Utyp ~ 1 — jma- Though formally (JTUJ) is not valid for the soft edge 
scaling, and should be replaced by an accurate formula for the mean resolvent 
involving second solution of the Airy equation Bi((), the above estimate 
works well for ( ^> 1). This gives for the mean value 

(c sc ) = = iW3 ( M _ A _^ 

where we have used C t ypC = y typ . We will indeed find below that (c sc ) ~ 
— \/C) is the most probable value of the curvature for ( ^> 1. In what follows 
we will set again c sc = c for brevity. 

Now we proceed to evaluating the full curvature distribution in the soft 
edge scaling limit. From (j!67p we have 

Su) a((, u) = dr cos (Vc + y ^ s w e- |w|r (49) 

/■oo / 3\ /-oo 3 

J dr sin f t( + _) e - |aj|r + J dr e r ^ 



+2 

Further defining 



r 3 . 



/?(c,C) = -< r r ete «««(C.w) (50) 
./-oo 27r 

and using (1491) we then find 

/3(c,C) = c 7 (c,C)-7 , (C,c) + 5(c,C) (51) 

where 

7(cC) = if cc^ + g^, y(C,c).| 7 (c,C) (52) 

and 

5(c,C)=^(-c)e- c ^. (53) 
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This yields the contribution to the level curvature distribution corresponding 
to the first term in 



V^(c,C) = [MQFM)-AmP{c,Q] (54) 



where /3'(C, c) = |/3(c, C). 

Note that relations (I151[ I152[) imply 



f°° 1 

/ rfc/3(c,C) = 7r^(C) and Ai(C)Bi'(() - Ai'(()Bi(() = - (55) 



which when used together with f[54|) ensure that the above piece contains the 
full normalization: dcV^\c,C) — 1- Now we further notice 

/%C)= / —e^\u\a(Cuj) (56) 



which implies from (1471) 

/°° Hi i f) f) 

— M$(C, u) = - [Ai'(()/3(c, - f ((, c) Ai(()\ = (c, C), 

(57) 

Similarly, 

°g M*(£ w) = ^ [/3(C, c) A<"(C) - Az'(C)/3'(c, C)] , (58) 
Taken together this gives second contribution to the curvature distribution: 

V {II \c,C) = ~ |/3(C,c) Ai"(() - ^(C)/3'(c,C) + P (/) (c,C)^lnp(C) 

(59) 

The derivative form ensures dcV^(c,Q = 0, as expected. Note also 

that using Ai"(() = CA(C) it is easy to check p'(C) = - [Ai(()] 2 . 
Note that from the definitions (I52II53I) we have: 



5'(c, C) = -c5(c, C), /(c, C) = c 2 7 (c, C) - Ai(() (60) 

which yields 

(3'(c,C) = -c(3(cX) + A l (C). (61) 
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Substituting this to (I5l"|) we get 

P«(c,C) = -/3(c,C) [cAz(C) + Ai'(0\ + A?{Q (62) 

which implies: 

§-V {I) (c,C) = ~^{PM [cAi(0 +Ai'(()}} (63) 
Similarly, we have 

/3(C, c) A»"(C)- Ai'(C)^(c, C) = /3(c, C) [c Ai'(() + Ai"(C)]-Ai(C)Aif(0 (64) 
resulting in 

V^(c,C) = -§- c {(3(cXncX)} (65) 
where we have denoted 

u(c, C) = cAi'(() + Ai"(() - [cAi(() + &'({)] \^ (66) 

The above formulae provide exact curvature distribution in the soft-edge 
limit. As they are quite complicated, it makes sense to work out several 
limiting cases of general interest explicitly. 

3.3 Ivestigating large curvature asymptotics: c — > +oo 

In this limit S(c,j) = 0. We can expand = 3- (l - ^ + ^ + ■ ■ J • 

Substituting to (|52|) we get 



= i E 4^ (2fc) (o = 4 ^(o + 4^"(o + 4^""(o + • • • (67) 

c c c c c 

fc=0 

which implies 

l\cX) = - 2 Ai\() + -^'\C) + ... (68) 

so that 

/3 7 (c, C) = c 7 (c, - Y(c, C) = -E V^^ (fc) (0 ( 69 ) 

fc=0 



16 



Now a simple calculation gives 

7f )(c,C) = -cP 7 (c,()Ai(0 - /3 7 (c,C)M0 + At 2 (() (70) 



= E ( -^r- [MO ^^(O - MO A<«(0] 

or explicitly, up to the terms of the order 0(c -4 ), using = ( Ai(() 

p«( c ,C) « ^ (MO 2 - C^ 2 (0) + ^' 2 (0 (7i) 
+^ (C [MO 2 - M0^"(0] - MO MO} 

Now, rewriting ( 166|) as 

Kc,C) = c (MO - M0^) + (MO - Ai '(0l^) (72) 
we get after straightforward manipulations 

^(cXHcX) = (MO - M0^) Ai(C) (73) 

{^ w (o (mo -^'(o^) - a-^co (^ao -mo^) 



fc=l 

and therefore 



^f ) (c,0 = -^{/3 7 (c,0(c,C)} 



| t ^{^- 1 H0(M0-M0^) 

-^ (fc ~ 2) (0 (^AO-MO^) 
«-|(M0 2 -CM0)-^ 2 (0 (74) 
+ ^ (^ 2 7§ ~ 2Cp(c) ~ 3A(C)A(C) ) + ■ ■ ■ 
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Adding up the two contributions we obtain the general expression 



^o=EV" n ' (75) 

k=2 

where 

V k = (2-k)M\0A^- 1 \0-M0A^ k \C) + (k-l)Ai'\0^ {k ' 2 \0 (76) 

-(* " ^l^j [Ai'(0Ai (k - 2 \0 ~ Ai(()Ai^\()} 

One finds that V\ — V2 — V3 — 0, so that the first non-vanishing term in 
the sum is V4, and therefore: 

r y (c, C) « i {-2Cp(0 + - 4A(C)A'(C)} + 0(c" 5 ), c -> +00 

(77) 

The above result confirms our intuition that the large curvature values occur 
when the two levels approach closely, hence the large-curvature tail exponent 
being dictated by the level repulsion mechanism is therefore universal, see 
[22]. For (3 = 2 this mechanism indeed predicts V(c — > 00) ~ c~ 4 , in full 
agreement with ( l77j) . 



3.3.1 Towards the bulk: ( — > —00 limit. 

In this limit we approach the bulk of the spectrum and it is natural to expect 
that the result will match the bulk curvature distribution (I40p . 
As C = — |C|j with |C| —> 00, we can rewrite fl52|) as 



2iri-co" c 2 + t 2 

* (78) 



2tt J_ x c2 + Kir' 2 

where we have changed the integration variable r — > |('| 1 / 2 r. The above 
integral is clearly amenable to evaluation by the saddle-point method. The 
s addle-point condition is r 2 = 1, hence r = ±1. Explicit calculation gives 
the leading-order contribution: 

7 (c, -C » 1) « Az(C) = 7 o(c, C) (79) 
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where we have used the asymptotic of the Airy function in the same limit: 

Ai(c)K 7^ C08 (^l 3/2 -l) < 80) 

Anticipating that we may need next to-the-leading order terms, let us con- 
sider the difference 7(0, () — 70 (c, (). Using the exact integral representation 
formula (j!5ip for Ai(() one can show that: 

< r\ in 1 ICI 3/2 r i\(\""(-^4) (t 2 - 1) i ,on 

1 -i r 1 j e «i- 



2t (c 2 + Id) 7-00 c 2 + Kk 2 
and performing integration by parts we arrive at the exact relation 

7(c» ~ 7o(c, C) = / 2 1 iAh / e V / 2 1 1 /" 1 2V2 dr 82 

2tt (c 2 + (CI) J-oo (c 2 + Kr 2 ) 2 

We again can evaluate the limit \(\ — > —00 in the above expression by the 
steepest descent method, and find 

7 (c, -C » 1) - 7o(c, -C » 1) « /U 2 | >,J A() = 7i(c, C) (83) 

(c 2 + Id) 3 

We conclude, that for —( ^> 1 we can use the following approximation 

7 (c, -C » 1) « - T -^ ]7r Ai(C) + 2 Ai'(C) + . . . (84) 
c 2 + |C| (c 2 + Kl) 3 

(We can actually check that assuming in f )84|) further c ^> |d and expanding 
reproduces the series ( |67|) .) Differentiating, we find to the same order 

7 '(C, -C » 1) « -^-r Az'(C) + 1 Ai(C) + TT^W ^"^ + • • • 

c 2 + |d (c 2 + ld) 2 {c+\Qr 

(85) 

Taking into account that S(c,() is exponentially small in such a regime, we 
neglect it and find correspondingly 
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(Ai'(()-Ai"(0 + cAi'(()) + ... 



(c 2 +|CI) 3 

Substituting this to dS2J) gives using p(() = Ai'(Q 2 - Ai"(()Ai((): 
Similar calculation yields 

£«« "< » « 53 + (?Ti(i? (cAm ~ Af(Q) 

and further substituting to ( 105]) gives after straightforward but lengthy alge- 
bra 

pM( C -C>l)«--^) C ^ (C) (89) 

[c > C>>ij c 2 + |C| (c 2 + |CI) 2 l " Jj 
+ (c2+ 1 |c|)2 {-2C (^ /2 (C) - C^ 2 (0) - aAi(C)MC) + §^« 2 (C)} 

Adding the two contributions gives the first nonvanishing term to be 

(90) 

Note that assuming c ^> |£| gives back exactly (ITT)) . In the present case 
— ( ^> 1 we, however, can further use that Ai(£) ~ ^1^174, ^4*'(C) 



P(c,~C»l) 



TCfAXC), p(C) « P'(C) ~ -^lCr 1/2 ~ ^y- We then see that 

the leading term is the first one, and arrive at the final expression: 



2 ttV(C) _ 2 if 

^(c 2 + 7T 2 p 2 (C)) 2 ~ [C 2 + K 2 ] 



y( C ,-C >!)«- — I Vo^ ^Z ro . o l2 » « = TP(0, (91) 



The formula (I9TT) precisely matches the bulk curvature distribution (HOI) , as 
was anticipated. 
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3.3.2 Away from the bulk: £ — > +00 limit. 

Again the idea is to apply the saddle-point method for ( — >■ 00. We start 
with the asymptotic analysis for the Airy function. We shall see that for our 
goal we need actually also next to the leading order corrections to the Airy 
function. Let us find them from the saddle-point method. We start with the 
representation valid for £ > 0: 

M0 = ^jy< 3/2 ^dT, C >0 (92) 

The saddle-points are r = ±i, and the relevant one is r = i as known from the 
asymptotic analysis of the Airy functions(see below). Denoting the steepest 
descent contour passing through r = i as T, we shift to it from the original 
contour running along the real axis. It is known that T runs asymptotically 
tangent to arg{r) = n/6 for K(r) — > 00, and arg(j) = 5n/6 for 9ft(r) — > —00. 
However, for finite 3?r (in particularly, close to r = i) we can consider the 
contour as running actually parallel to the real axis and parametrise V as 

T = i + (^374- Tms S ives 



1 _lr3/ 2 T/ \ t / \ f°° e -4+^ dv 



W^C 1 J-00 V27T 

(93) 

The corresponding leading-order asymptotic for the Airy functions is ob- 
tained by replacing J (a) = 1: 

^(C)^^o(C) = ^ T7I e-K 3/2 , (94) 

and a similar calculation for Bi gives 

Bi(C) « Bi (C) = ^T7I^ 3/2 , C»l- (95) 

These relations imply to the leading order: 

Ai'(C) « -VCAi Q (0, Bi'(() w y/CBi (C), ( » 1 (96) 



which indeed gives the correct Wronskian, cf. fl55|) . However, for our present 
goals we will need to account for subleading terms as well. To see this 
fact we can consider the mean density p(Q = Ai'(Q 2 — (Ai 2 (() featuring 
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in the calculation. Using 



immediately shows that the leading order 



vanishes, and we indeed need to go to the subleading order by expanding 
J (a^> 1) = 1 — 4r + 0(a~ 4 ). This gives for the Airy function: 



Ai(c)«A<o(o(i-^+o(r 3 : 



(97) 



After differentiation, using (J94 



^'(C)^-Az (C)C 1/2 i 



11 1 

16 

With this expressions we find using p'(C) = ~Ai 2 (() that 



+ 0(C 3 ) , C»l (98) 



p(0 



1 [Az (o] 2 (i + o(r 3/2 )), ^«- > /c^i + 



2v/C 



MO 



16C 3/2 



(99) 



with some yet unspecified coefficient b. Consider now the combination (166 
conveniently re-arranged 



i/(c,C) = c 
We have 



+ CA(C)-^(C)ig 



(100) 



CAz(C) - ^'(0^4^ ~ ^o(C)-^ + 4) 



M0 



16VC 



so that 



u(c, C) « Ai (C)— [c(b - 4) - v^(6 + 4) 



(101) 



Therefore, we have asymptotically to the leading order 

7^)(c, C) « -Aio(0^§- c {/% C) [c(6 - 4) - ^C(6 + 4)] } (102) 

and to the same order 

pW(c,C) = -/3(c,C) Mi(C) + Aif(0) + A?{Q 
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« 0(c,()Ai o (O(V(-c) +Aig(C) (103) 

It remains to find asymptotic for /3(c, £). To this end, we perform a similar 
analysis for j(c, £). We write 

^0 = ^]_Jre y a) ^^ (104) 



Attic 



r - ic/vT r + z'c/ vT. 

For |c| < deforming contour from the real line to the steepest descent 
contour V incurs the contribution from the pole at r = i\c\/y/(, whereas for 
|c| > yf( the pole is above T. We therefore have the exact identity: 

7 (c, C > 0) = 7p (c, C > 0) + 7 r(c, C > 0) (105) 

where we denoted the pole contribution 

7p (c, C > 0) = ±-^ + ^6(^/C - |c|) (106) 

and 

Performing the integral by the steepest descent method gives ( using ( 1941) ) 

7r (c,C > 0) « Aio(C), C » 1, |c| ^ (108) 

c - C 

It turns out that for |c| < ^/C t ne term 7 p gives much larger contribution 
than 7r- To verify this we rescale c = CyfC, and consider c to be of the order 
of unity. For 7 p we arrive at the expression: 

7 P (c, C > 0) = _^ e -C' a (/(M^(l _ | g |), /(„) = „ - in 3 (109) 

As /'(u) = 1 - u 2 > for u e [0, 1) we see that f(u) < f{l) = 2/3, hence 
7 P (c, ( > 0) ^> ^TT/^e - ^ 372 . At the same time using the asymptotics of the 

Airy function 7r ~ 2 ^J^ 5;/4 e~i^ 3/2 < ^r^e - ^ 372 <C 7 P , proving the statement. 
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Thus for |c| < \/C we can approximate j(c, £) ~ 7p(c, C)> which upon 
differentiation gives 

7;(c,C>0) = -V^l c ^(v^-|c|) 

and further 

c lp (c, C > 0) - 7 ;(c, C > 0) « e^ c+ ^ c3 0(c)0( - c) 
and finally adding (|53|) we obtain 

/% C > 0) « /3 p (c, C > 0) = c 7p (c, C > 0)- T ;(c, C > 0)+5(c, C) - f 

(110) 

which is the leading approximation ( exponentially dominant) in the whole 
domain c G (— oo, a/C)- 

It is useful to check that such a precision is sufficient for the relation ( 155]) 
to hold: 



/3(c, C) c?c 



/V, 
e^+^cfc^ 7rSi(C) 
■oo 



111^ 



which is checked by performing the integral (after change c — > a/C c ) by the 
steepest descent method, and using asymptotic ( 1951) for Bi((). 
Substituting (lllOp to (I103P gives the leading order expression: 



7f )(c,C) « Aio(C) e- Cc+ 3 c3 (v/C - c), c G (-oo, vf) 
Differentiation in the domain c G (— oo, y^) gives 

^ p (c,C) = (c 2 -C)/3 P (c,C), cG (-oo,VC) 
which finally shows that to the leading order 

VW(c, C) « Aio(C) ^-^(C-c 2 )^ [c(6 - 4) - VC(& + 4) 



;ii2) 



113) 



c G (— oo. 



(114) 

Adding together we get the full curvature p.d.f. for £ 1, and c G (— oo, y/C): 



(Vc-c){i + i(i + ^ 



—(6-4) -(6 + 4) 



(115) 
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To understand the structure of the above expression notice that 
Ai (() e^ c+ ^ = —l—e^-^-l^ 2 = _l_ e 4(VC- C )( C +v/C) 2 (116) 

so that the curvature p.d.f. has a very sharp maximum (with the height 
£ 1//4 and widths C" 1 ^ 4 ) around c = — a/C Introduce correspondingly the new 
random variable x = (c + v / C)v / 2C 1/ ' 4 ) so th & t c = ^ 1/4 — y^. 

Ai (C) e" Cc+ ^ 3 = ^T7I e "^ + ^^ 

Remembering that the p.d.f. of x acquires the extra Jacobian factor \/2C 1//4 , 
we see that the probability density 

7>,(,,C » 1) « ^e-*^* (l - ^) {l + O (^)} (118, 
tends to the standard gaussian distribution -75=e 2 for fixed x and C — > 00. 

u V 2-7T 

Now it remains to consider the case c > -\/C- Then the only contribution 
comes from the term 7r(C)- Using precisely the same method as in £ — > —00 
limit we can arrive at the full analogue of the exact identity (1821) : 



7r(c , _ _1_ AiK) = £ e -c-(-=f)_L_ rfT (119 ) 

and performing the integral by steepest descent further see 

7r(c, C) « Ai(C) + Ai(C) + . . . (120) 

Noticing that |£| = — C for £ < 0, we see that (11201) is precisely the same as 
hence the curvature p.d.f. will be again given by analogue of 



Vie > VI » 1) « {-2CP(0 + ~A 2 ^ - 4A(CK(C)} + • • • 

(121) 

Using the asymptotic formulae ( 1971) . (1991) we find that the leading term in 
(11211) cancels, and the result is of the order of Ail(Q(^ 3 ^ 2 ~ or explicitly 

. rz n 2.B P\C) B 1 4^3/2 

V(c>^>>l)~ w - w /-f* w - w — 2 e-*<, (122) 
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where the constant B of the order of unity is left undetermined. 

The crossover between the two regimes c < a/C an d c > a/C happens over 
the domain |c — y/£\ ~ ^74. Indeed, take c = a/C + > where /3 is of the 
order of unity. Substituting this for /3 > into 7r from (j!08p and using 
asymptotic for Airy function at C ^> 1 we see that 

1 _2 f 3/2 1 _2a3/2 /ino\ 

■<^ 2vc(c-cw< e ' K v^ es (123) 

On the other hand, the pole contribution f |106[) for c ~ a/C approaches the 
value 

>^e-fC- (12 4) 

that is the two contributions are of the same order. It would be interesting 
to find the exact crossover expression for arbitrary fixed /3, and we leave it 
for further investigation. 
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Appendix B: Bulk asymptotics of Hermite polynomials, Cauchy trans- 
forms, and the kernels involved. 

As is well-known ( see e.g. [T7]) the Hermite polynomials Pk(x) have the 
following integral representation: 



( \ N 

V Z7T 



(-i) N +™I N+n (x)+i» +n I N+n (x) 



N+n - 



where 

lN+n{x) 



dq q 



N+n e -^{q-ix) 2 



N — 

e 2 



N+n e -^-q 2 +iNxq 



(125) 



(126) 



For real \x\ < 2 we can use parametrization denoting x = 2cos0 and by the 
steepest descent method (s.p. at q — ie~ l< ^) find p2] the following large— A 
asymptotic behaviour: 



lN+n{X) 



■N+n 



71 



Asin< 



3 f cos 2^ -i{n+\)<j)+i^-iN9{(l>) 



i-- sin (20) 
(127) 
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so that 

Piv(x) f« \ -r—r e 2 



sm i 



cos 2 * cos - ^ + iV0(0)^ (128) 



Differentiating over 2 , picking up the leading terms proportional to N, and 

d<t> 



usin s t = -2ib and ^ = 2 sin2 (0) we s et 



ax 2 sm y sm 

- sin(20) cos ( ^0 - - + iV0(0) ) - 2 sin 2 (0) sin ( ~0 - - + N0(<j>) 



= Np N _ 1 (x) (129) 

which in fact is easy to show to be the exact relation for Hermite polynomials. 
Similarly 

Pn-i(x) w *Eef cos2 ^cosf~0-^ + iW(0)Y ^p N ^(x) « jV^_ 2 ( 
V sm0 \ 2 4 J ax 

(130) 

This shows that the kernel W^i(x, x) from ( 125]) can be written as 

Wifox) = ATfpAr-i^iv-i^) -^(^-2^)] ~ 2iVsin0e Arcos2 * (131) 

where we have used the identity cos 2 (A) — cos (A + 0) cos (A — 0) = sin 2 
for any A. This implies 

W 2 (x,x) = —Wx(x,x) « (-2iVsin (20))( - Wi(ar,z) = 2A^cos0iyi(x, 

ax 2 sin 

(132) 

Finally, recalling the mean eigenvalue density p(x) = ^V^-? we re- 
cover (1371). 



Now, denote the Cauchy transforms of the above Pk(x) as h k (x) = fk(x). 
As is shown in [T7] for real x and e holds the exact integral representation 

f N+n {x + ^e) = -{^^^^J^e-^^^^dq (133) 
where we denoted s e = sign(e). Comparing with (1126ft we conclude: 

f N+n (x + it) = -(-i) N+n s? +™- 1 ,/^ e -f ^ +i ^I N+n (xs e +i\e\) (134) 

V Z7T 
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In our applications we will need e = i^y (see eqs. (??) and (?? |23ll24l) ) 
when considering the large- N asymptotics. To this end, the formula (1 127ft 
obviously implies in the large— N limit to the leading order, with real ( = 
0(1) 



(135) 



where 



Sx 



i c 



6' (0) = 1 - cos (20) = 2sin 2 0. 



dx 2 sin <p N ' 

Taking into account sin0 = | a/4 — x 2 = 7ip(x) we rewrite the above as 



I N+n (x + jj) « ijv +n (x) f 



£(cos <j)+i sin < 



Jjv + n(z) e' 



C(f+mp(x)) 



(136) 



which retains its validity when we replace real ( with purely imaginary ( 
\ y as is needed in our application (remember ?/ > by definition). 



(137) 



Combining this with (I134p gives needed leading-order asymptotic of the 
Cauchy-transform: 



f N+n (x + i^-y 



2tt 



(138) 

Now we can calculate the asymptotic of the kernel F\ (x, x + ij^y) from (123 p . 
Actually, it is more convenient to consider a slightly more general case: 

Fl ( x + ~k' x + = f N { x + Pn - 1 { x + jj) ( x + i-^y) p n ( x + /; 



(139) 



N 



_^_^N-l s N-2 e -fx 2 e -icu^y-np(x)\u\y 

Assume first u > 0, then 



-lS u I N (xS u )p N -i (x + JjJ - In-i(xS u )p N (x + ^ 



-is u I N (xs u )p N -i (x + j-j - In-i{xs u )p n (x + j- j 



-U N (x)p N -i (x + -^J - I N -x(x)p N (x + j- 
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or, using (j!25p . (I136p and the relation I N (—x) = In(x) for any real x, we 
have (up to the overall factor y|^) 

1-U N (X) {(-if-'lN^Xy^ + C.C.) - I N ^{X) {{-i) N I N {xy^ + C.C)} 



oc e 2 



-i N elV+iKp^ri 



[I N (x)I N ^x(-x) + I N (-x)I N -i(x)] 



Now assume u < 0, then similar calculation gives: 



is u I N (xs u )p N -i (x + %=) - In-i(xs u )p n [x+ ^- t 



N 



N 



N 



H N (-x)p N -i (x + — ) - I N -!(-x)p N [x + — 



N 



N 



= _ ( _,)^ e f™W^_ [i N { x ) Ilf _^- x ) + I N (- x )I N _x(x)} 
Thus we can conclude that 

-is u I N (xs u )p N -i (x + j^j - I n ^i{xSu)pn [X + j- J 



(140) 



I N {x)I N -x(x) + I N (x)I N -i(x) 
which gives for the leading-order kernel asymptotics: 



Fx [x + —,x + t—y 



N 
2tt 



~ £jtl_ g~T :r2 ef (v-iuy)-np(x)\u\y-iTrp(x)ri 



I N (x)I N -i(x) + I N (x)I N -i(x) 



(141) 

Using (JI27D we find I N (x)I N _x(x) = ijr^e Ncos ^- l(t> , so that finally we 
find for the kernel: 



CO 



F 1 [x,x + i—y 



N„2 



x —Uk)%y—Trp{x)\w\y N cos 2<f> 



I e 2 e 2 



ie e 



iV„-iu>§ y-irp(x)\u\y 



(142) 

equivalent to the equation (1331) . 

Finally, the kernel F2 (2, x + ij^y) is obtained by differentiating ( 1141H 
over r), setting 77 = 0, and multiplying by factor N. which yields the relation 



F 2 [x,x + i—y 



N 



~x , ."1 / 

- - «s w 7rp(x)J F x \ x, x + 1— 1 



(143) 



which is equivalent to ( 136)) . 
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Appendix B: "Soft edge" asymptotics of Hermit e polynomials, Cauchy 
transforms, and the kernels involved. 

We start again with 

I N+n (x) = e N ^J N+n (x), J N+n (x) = dqq N+n e-^ 2+m ^ (144) 

Jo 

and replace the contour with the sum of two contours [0, i] i + oo], so 
that, correspondingly, J N+n (x) = JjP +n (x) + J^P n (x). In the first contour 
we parametrize q = ip, p G [0, 1], so that 

J^ix) = i N+n+1 I dpp N+n e% p2 - Nxp (145) 
Jo 

and in the second contour we put q = i + t,Wt > 0, so that 

POO 

jM n ( X ) = dt{i + t) N+n e -^(t + i? + iN x{ t + i) (M6) 

Jo 

We will be interested in the regime x = 2 + ^§73 j with ( of the order of unity, 
and N ^> 1. Let us start with rewriting (j!46p as 

roc 1 j.2 

jW n {x) = e~l N / dt ( i+ t) N+n e^^e- NC ^, C(t) = ^-it-ln (i + t) 
Jo 2 

(147) 

The saddle-point equation is ^ = = t — i — which has the only solution 

t = 0. Expanding for £ <C 1 gives C(t) fa — lnz — iy + 0(t 4 ). Introducing 
the scaled variable: t = we easily find to the leading order 

jN+n 1 poo 3 

J^s e 2Ne ~ N ^(0, a dO = l dre^V (148) 

Similarly, in the integral (j!45p we make the substitution p — 1 — ^773, with 
r G [0, iV 1 / 3 00] and after replacing x = 2 + -^73 and expanding for 
r <C iV 1 / 3 find to the leading order 

■N+n+l J roc 3 

jgU^wi^-e-i^e-^faaCC), « 2 (C) = jf dre^V (149) 
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Combining, we finally have for x = 2 + 



C 

AF273 



■N+n x 

lN + n(x)^j^e* N e N ^a((), o(C) = a x (C) + 102(C), V( G C (150) 
Note: Introduce the two functions: 

Ai(C) = ^ drcos^rC + y) (151) 

and 

52(C) = - / <*r sin rC + — + - / dr e T< ~ (152) 
n Jo V 3/ ttJo 

which are the two standard linear independent solutions of the Airy equation 
/"(O - C /(C) = 0. Obviously, o(C) = vr [Ai(() + iBi(()}. 

After substitution to (1125)1 this yields (for real C and Wn <C N) to the 
leading order : 

PN + n (2 + j^j « v/27iV 1/6 e*" e^(C), (153) 

We see that to this order there is no dependence on n, which will result 
in vanishing of the corresponding kernel. To find beyond-the-leading order 
corrections we will use the exact recursion: p^^i(x) = jj-^Pn(x). In the 
"soft edge" scaling regime we, correspondingly, have 

P»-i ( 2 + ivWa) = WHTC PN ( 2 + Wi) ' (154) 

Pn - 2 ( 2 + tU) = 4^ie PN ( 2 + -m) (155) 



which results in 



P»i 2+ Jjm) *A N e Ni <Ai(C), A N = N 1 / 6 e 1 * 



N 



PN-2 



PN-i (2 + J ^j^A N e^ 
( 2 + ^H^ 



mo + -j^m) 



(156) 



mo + j^*m + j^-s^"(o 
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Remembering the relation (11311) we have for the kernel W\{x,x) at x 
2 + 



AT2/3 



Wi(x,x) = N [p 2 N ^(x) - Vn {x)vn-2{x)] « 2 7 riV 2 / 3 eV 7V3C p(C) (157) 

where we denoted p(() = Ai'(() 2 — Ai(()Ai"(() (which indeed is proportional 
to the mean GUE eigenvalue density in the "soft edge" regime). Then in view 
of the exact relation W 2 (x,x) = ^Wi(x, x) = N 2 ^ 3 -^Wi(x,x) it follows to 

the leading order in N 3> 1 W2(x,x) ~ 2NWi(x,x) for x = 2 + 



TV 2 / 3 ' 



This 



is indeed compatible with the "bulk" relation fl3T|) W 2 (x,x) ~ NxW\{x, x) 
as now x ~ 2. We shall see however that for our goals we need to keep the 
corrections to this approximation: 



W 2 (x,x)^2nN 2 / 3 e N e 2N ^ 



2Np(0 + N 2 ^±p(0 



(158) 



so that 



W 2 (x,x) _ r o /o 1 d ... 

N + iV 2 / 3 ^^ — p(C), a; = 2 + 



2W 1 (x,x) "' ' " 2p(C)c/C rV ^' ~ " ' ^ 2/3 (159) 
In what follows we will also need a similar recursion for In{x), which is 
simply /jv_i(x) = — -^-^In(x). It can be easily derived from the integral rep- 
resentation (11261) by using the identity ^ ^e~^ ( - q ~ lx ^ j = — i-^ (e 
and integrating by parts. This implies 



in 



l N-l 



c 



2 + 



AT2/3 

c 



A N e N ^a{(), A N 



;N 



C 2 



JV 



'7V-2 



2 + 



AT2/3 



AT2/3 
«(-/) 2 i. Y r- V ^ 



fl (0 + ^'(0 



(160) 



«(0 + ^a'(C) + ^a"(C) 



Now, we use (II 34ft for x = 2 + -^73 and e 



jyi/3 

my , 
7V273 ■ 



AT+n 



2 + 



iV 2 / 3 



-0 



iV+n AT+n-1 




x e 



jv2/3 ; 7 



7V+n 



2s w + 



(s u + i\u\y 

jV 2 /3 



(161) 
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Using fll50p we get the leading order expressions for uj > 0: 

fs (2 + - B N e- N ^v) a((+iuy), B N = e~l N u>0 

(162) 

where 

pOO 2 pOO 2 

a(( + iuy)= dre !rC+,: Te"^ T + « / dr e T ^^ e^ T (163) 
Jo Jo 

Similar calculations can be performed for u < (this requires to know 
asymptotics of In(x) around x = —2) with the result 



f N (2 + « e-^*«-*«) 5(C+^), = e"^ u, > 

(164) 



where 

a(C + iuy) 



poo 3 /*OG 3 

/ dr e- iT ^ iT - e^ T - i / dre^^e^ (165) 
Jo Jo 



Combining these two formulas for any real u we can write 



A' 



^2/3 J "vs,~v, -iv 

(166) 

where 



Jo Jo 



iV-1 



2 + 



AT2/3 



Q t (C,w) + -^ i7 3a , (C,w) , (168) 



where the dash stands for the derivative over £. 
Now we substitute (II 68ft and (1156H to the kernel 



Fi ( x, x + 



iuy_ 

AT2,3 



fN\X + 



iuiy \ 



p N -i(x)-f N . 



x + 



iury_ 

jV 2 /3 



Pn{x) 



x = 2+ 



AT2/3 



(169) 
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and obtain 

m (2 + -^- 2+ C + lUJ y ) a 
e -iv e -i^ u ^ [A . (C) a , (c w) _ A , (c) a(c u)] ^ (17Q) 
A similar calculation gives also: 

I [Ai(C) «'(C, w) - MO «(C, w)] + ^73 [Az"(C) a(C, w) - a'(£ w) MO ] } , 

(171) 

We see that to the leading order F 2 (x, x + -0^) ~ -/VFi (x, x + j^tfs) in the 
regime x = 2 + -^73 , in full agreement with the x — > 2 limit of (136]) . We shall 
see however that for our goal we need the full expression ( 11711) . 
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